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ABSTRACT 

We present the first dynamical determination of the binary parameters of an eclipsing 
SW Sextantis star in the 3-4 hour orbital period range during a low state. We obtained 
time-resolved optical spectroscopy and photometry of HS 0220-1-0603 during its 2004- 
2005 low brightness state, as revealed in the combined SMARTS, IAC80 and Ml Group 
long-term optical light curve. 

The optical spectra taken during primary eclipse reveal a secondary star spectral 
type of M5.5±0.5 as derived from molecular band-head indices. The spectra also pro¬ 
vide the first detection of a DAB white dwarf in a cataclysmic variable. By modelling 
its optical spectrum we estimate a white dwarf temperature of 30 000 ± 5 000 K. 

By combining the results of modelling the white dwarf eclipse from ULTRA¬ 
CAM light curves with those obtained by simultaneously fitting the emission- and 
absorption-line radial velocity curves and /-band ellipsoidal light curves, we measure 
the stellar masses to be Mi = 0.87 ± 0.09 Mq and M 2 = 0.47 ± 0.05 Mq for the white 
dwarf and the M dwarf, respectively, and an inclination of the orbital plane of * « 79°. 
A radius of 0.0103 ± 0.0007 R© is obtained for the white dwarf. The secondary star in 
HS 02204-0603 is likely too cool and undersized for its mass. 

Key words: binaries: close - stars: individual: HS 02204-0603 - stars: fundamental 
parameters - novae, cataclysmic variables 


1 INTRODUCTION 

Cataclysmic variables (CVs) are binary systems that con¬ 
tain late-type (K-M) companion stars that fill their Roche 
lobes and therefore transfer mass to their white dwarf (WD) 
primary stars. Nova-like (NL) variables are CVs for which a 
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nova or dwarf nova outburst has never been observed. The 
orbital period distribution of CVs shows a signihcant paucity 
of systems in the ~ 2-3 hour range, the so-called period gap 
(see e.g. [Gansicke et al.||2009] and references therein). CVs 
with orbital periods above the period gap are expected to 
have M-dwarf secondary stars up to a period of about ~ 5.8 
h, while secondary stars of earlier type are to be found for 
longer orbital periods (Knigge 20061. It is also just above 
the period gap, between 3 and 4 hours, where a large popu- 
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lation of CVs with extreme behaviour is found (Thorstensen 


et al. 1991 Rodriguez-Gil et al. 20071. There is mounting 


evidence that the NLs that populate this orbital period in¬ 
terval harbour intrinsically very bright accretion discs and 
very hot WDs (Townsley fc Bildsten|2003 Araujo-Betancor 


et al.|2005 Townsley fc Gansicke||20TO I that current evolu¬ 

tionary theories are unable to account for. In addition, at 
least 50 per cent of them belong to the SW Sextantis class 


I RodrIguez-Gil et al. 20071. Their large intrinsic accretion 
luminosities completely outshine both the WD and its com¬ 
panion most of the time, making direct observations and 


accurate binary parameter measurements impossible (Cia- 
rdi et al.|1998 1. However, these CVs are occasionally caught 


fading towards states of greatly diminished brightness, or 
‘low states’. During the low states, these CVs are 3-5 mag 
fainter than their normal high state. They can stay at this 
level for days, months or even years before returning to the 
high state. Low states appear to be a feature of several types 
of CV and, in addition, they seem to occur independently of 
the magnetic field of the WD: they are observed in practi¬ 
cally all known strongly magnetic systems without discs (po- 
lars), some intermediate polars, and many weakly-magnetic 
CVs, including a large fraction of NLs and a number of dwarf 


novae (e.g. Schreiber, Gansicke & Mattei 2002 Manser & 
Gansicke||2C)14 1. 


Although there is broad agreement that low states are 
the consequence of a reduction in the mass transfer rate 
from the donor star, the exact mechanism is unclear fLivio| 


& Pringle 119941 and King & Cannizzo (19981 proposed 


an accumulation of large starspots close to the inner La- 
grangian point (LI) as a way of inhibiting Roche-lobe over¬ 
flow, but Howell et al. (20001 put the case of the polar ST 


LMi forward to challenge this scenario and proposed that it 
is changes in the level of magnetic activity of the secondary 
star that may be behind the observed high-/low-state tran¬ 
sitions. This idea had been already explored by |Bianchini| 
(19921, who invoked solar-like magnetic cycles with time 


scales of about a decade. 

In the low state, only gas provided by the magnetic ac¬ 
tivity or stellar wind of the donor star would then be avail¬ 
able for accretion (e.g. [Hessman, Gansicke &: Mattei||2000p . 
In this respect, entanglement of the magnetic fields of both 
stars has been proposed to explain the line emission patterns 
observed in polars during low states (e.g. Kafka et al.|2008 
Mason et al.|2068 l. 


The study of the GV population with orbital periods 
between 3 and 4 hours is of essential interest to the develop¬ 
ment of the theory of CV evolution. Particularly, those NL 
systems closer to the upper boundary of the period gap at 
an orbital period of about 3.18 hours (a limit set observa- 
tionally by Knigge|2006 1. According to the standard stellar 
structure codes (e.g. |Chabrier fc Baraffe|1997| ), the M-dwarf 
secondary stars in these CVs are expected to undergo a ma¬ 
jor change in their internal structure, developing into a to¬ 
tally convective configuration at a mass of ~ 0.35 Mq. The 
exact value of this boundary depends on the choice of sev¬ 
eral parameters, such as the metallicity, opacities and mix¬ 
ing length (e.g. [Chabrier fc Baraffe||1997[). We can find in 
the literature values ranging from 0.27 ( Gabriel|1969 ( to 0.4 
Mq (Cox fc Giuli|[l9^ I. For low-mass stars Reiners & Mo- 


|hanty| ( |2012[ ) pl^e 0.3 Mq as the fully-convective limit. In 
the case of detached white dwarf/main sequence (WDMS) 


binaries a drop in the fraction of post-common envelope bi¬ 
naries (PCEBs) relative to wide binaries is observed at a 
mass of ^ 0.25 Mq, which is an independent observation in 
support of the disrupted magnetic braking mechanism in¬ 
voked to reproduce the observed orbital period gap in CVs 


I Schreiber et al. 20101. On the other hand, in the population 


synthesis study of CVs by Howell et al. (20011 their zero-age 


main sequence models have a fully convective configuration 
at a mass of ~ 0.34 Mq, while the semi-empirical donor star 
sequence for CVs of Knigge (20061 suggests a (fully convec¬ 
tive) donor mass of 0.2 Mq at the upper edge of the period 
gap. These authors warn against the lack of fundamental 
input in the form of accurate masses, especially in the 3-4 
h orbital period range. In fact, the donor masses used in 
this period interval come from indirect methods such as the 
mass ratio (q)-superhump period excess (e) relation refined 
by Patterson et al. (20051, and the a priori assumption of a 


constant WD mass of 0.75 Mq. Moreover, the system used 
to refine the fit to the q-e data, thought to have the largest 
value of the superhump excess (e = 0.094), is BB Doradus, 


(Rodrfguez-Gil et al. 2012 Schmidtobreick et al. 20121 for 


which an accurate determination of its stellar masses is un¬ 
available. 

Another essential observational input to the theory is 
the time-averaged mass transfer rate, that has to be com- 
pared with the predicted secular mass transfer rates (see 
Knigge, Baraffe fc Patterson|20lT I. A medium-term (10^ — 
10° yr) mass transfer rate may be derived from measure¬ 
ments of the effective temperature and mass of the accret¬ 
ing WD ( [Townsley &: Bildsten|[2003 Townsley & Gansicke 
20091, which again emphasises the need for accurate binary 


parameters. 

In brief, no direct measurement of the stellar masses in 
any SW Sex star in the 3-4 h orbital period range exists thus 
far. However, the CV evolution theory demands appropriate 
input in the form of dynamical mass solutions and other bi¬ 
nary parameters. This is especially true for these CVs close 
to the upper boundary of the period gap, in which the donor 
stars may be about to experience the predicted major struc¬ 
tural change in their stellar interiors. 

In this paper we present the first attempt to measure the 
fundamental parameters of a SW Sex star in the 3-4 h period 


range, the eclipsing system HS0220-1-0603 (Rodrfguez-Gil 


et al. 20071, with no initial assumptions about the binary 


system. 


2 LONG TERM LIGHT GURVE 

The long-term V-band light curve of HS 0220-1-0603 is pre¬ 
sented in Fig. The photometry spans 9.92 years and 
comprises data provided by Spanish amateur astronomers 
of the Ml Group, the 1.3-m SMARTS telescope on Cerro 
Tololo, Chile (operated by the SMARTS Consortium), and 
the 0.82-m IAC80 telescope at the Observatorio del Teide 
on Tenerife, Spain. The apparent V-band magnitudes were 
computed relative to the comparison star USNO-A2.0 0900- 
00554112. Conversion of its Ba 2 .o and Ra 2 .o magnitudes into 
the Landolt standard systerrQ resulted in a V-band magni- 

^ See http://www.pas.rochester.edu/-emcmiajek/memo_ 

USN0-A2.0.html 
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Figure 1. Top panel: Long-term, V-band light curve of 
HS 0220+0603 showing 9.92 years of photometry from 2003 Jan 
26 to 2012 Dec 26. Data from the Spanish amateur astronomers 
of the Ml Group (filled circles), the 1.3-m SMARTS telescope on 
Cerro Tololo (open circles) and the 0.82-m IAC80 telescope on 
Tenerife (open squares; first two data points) are shown. Data 
taken during eclipse are also included in the plot. Bottom panel: 
Zoomed version of the top panel showing the low state data only. 
The vertical dashed lines mark the dates of the spectroscopic, 
photometric and XMM UV observations detailed in Section]^ 


tude of 12.3. The nearby star USNO-A2.0 0900-00553941 
{V = 16.5) was used as the check star. 

The Ml Group reported HS 0220+0603 entering a low 
state on 2004 September 19. The system was a,t V — 16.8, 
nearly a magnitude fainter than the previous IAC80 tele¬ 
scope measurements. Almost a month later, on 2004 Octo¬ 
ber 24, the system settled in the low state at H = 19.5 after 
fading at a rate of ~ 0.08 mag/day. The egress from the 
low state occurred shortly before 2005 July 19, when the 
system was observed at H = 19.2. It then brightened at a 
relatively slow rate of ~ —0.05 mag/day and reached its high 
state brightness level at 1/ « 16 around 2005 September 18. 
This places a lower limit to the low state length of about 
1 year. However, just after 2005 September 18 it started to 
fade again at a rate of ~ 0.03 mag/day. HS 0220+0603 re¬ 
mained at V ~ 16.7 for about 64 days before starting its 
final rise to the high state at the same rate of ~ —0.03 
mag/day. All the ingress and egress slopes reported here 
have to be taken with caution due to the scarce sampling 
of the long-term light curve, but the occurrence of the sec¬ 
ond fading looks convincing. Further, during this secondary 
plateau the system showed occasional brightenings with an 
amplitude of ~ 0.5 mag, similar to those seen in BB Dor in 
the low state (Rodn'guez-Gil et al.||2012|. Finally, the post- 


low state monitoring reveals a steady brightness decline at 
a rate of ~ 3 x 10“^ mag/day. 


3 DATA ACQUISITION AND REDUCTION 

3.1 Optical spectroscopy 

3.1.1 Very Large Telescope 

We obtained four 600-s spectra with the ESO 8.2-m Very 
Large Telescope (VLT) at Cerro Paranal on 2004 Novem¬ 
ber 20, shortly after the onset of minimum light with the 
FORS2 spectrograph on the Antu (UTl) telescope. For 
the observations the spectrograph was equipped with the 
600RI grism, the GG435 order sorter filter and the MIT/LL 
2048 X 4096 pixel CCD detector. With a 0.7-arcsec slit, 
this delivered a spectral resolution of 3.5 A (full-width at 
half-maximum, FWHM) over the 5015 — 8325 A wavelength 
range. For flux calibration we used spectra of the flux stan¬ 
dard LTT 377. 


3.1.2 William Herschel Telescope 


We carried out time-resolved spectroscopy of HS 0220+0603 
on 2005 January 3 with the 4.2-m William Herschel Tele¬ 
scope (WHT) at the Roque de los Muchachos Observatory 
on La Palma. The double-armed Intermediate dispersion 
Spectrograph and Imaging System (ISIS) was fitted with the 
R600B and the R316R gratings, and the 2048 x 4200 pixel 
EEV12 and the 2148 x 4700 pixel Marconi CCDs in the blue 
and the red arm, respectively. This setup and a 1.2-arcsec 
slit width (~ 1.5-arcsec seeing conditions) granted access to 
the wavelength intervals 3700 — 5010 A and 6110 — 8930 A 
at 1.8 and 3.3 A resolution (FWHM) in the blue and red, 
respectively. The exposure time was 1200 s and we regularly 
took spectra of arc lamps in order to account for instrument 
flexure. A flat field image was obtained after every two sci¬ 
ence exposures of the target in order to correct for fringing 
in the red arm CCD. Flux calibration was performed using 
spectra of the flux standard Feige 34 taken with the same 
setup. We used the same spectrophotometric standard star 
to remove the telluric features. See Table [T] for details of the 
observations. 

We reduced both the VLT and WHT spectroscopic 
data using standard procedures in IRAI0 After subtract¬ 
ing the bias level, the images were divided by an average 
flat field that was normalised by fitting Chebyshev func¬ 
tions of high order to remove the detector specific spec¬ 
tral response. We then subtracted the sky contribution and 
optimally-extracted the final spectra following the method 


described by Horne (19861. For wavelength calibration, we 
fitted a low-order polynomial to the pixel-wavelength arc 
data. In calibrating the WHT data, we obtained the wave¬ 
length solution for each target spectrum by interpolating 
between the two nearest arc spectra. All subsequent analy¬ 
sis was performed with the MOLL’ij^package. All the spectra 


^ IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories. 

® http://deneb.astro.Warwick.ac.uk/phsaap/software/ 
molly/html/INDEX.html 
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Table 1. Log of optical observations. 


UT Date Coverage 

Filter/Grating 

Exp. 

Frames 

(h) 


(s) 


VLT, FORS2 spectroscopy 



2004 Nov 20 0.54 

GR1S 600RI 

600 

4 

WHT, ISIS spectroscopy 



2005 Jan 03 3.92 

R600B/R316R 

1200/1200 

12/12 

NOT, ALFOSC photometry 



2004 Dec 04 3.62 

7#11 

40, 60 

181 

2005 Feb 02 2.22 

7#11 

30, 40, 60 

107 

2005 Feb 05 2.35 

7#11 

40, 80 

142 

WHT, ULTRACAM photometry 



2005 Aug 09 4.27 

u', g', i' 

4.95 

3089 

2005 Aug 10 3.93 

i' 

4.95 

2841 


were corrected for interstellar reddening using a colour ex¬ 
cess of E{B — V) = 0.047, corresponding to Av = 0.145 mag 
(from the extinction map of Schlafly fc Finkbeiner|[20TT 1. 


3.2 Optical photometry 

3.2.1 NOT I-band light curves 

We used the Andaluci'a Faint Object Spectrograph and 
Camera (ALFOSC) equipped with the 2048 x 2048 pixel 
EEV chip (CCD #8) on the 2.56-m Nordic Optical Tele¬ 
scope (NOT) on La Palma to obtain time-resolved, /-band 
photometry of HS 0220-1-0603 on 2004 December 4 and 2005 
February 2 and 5. The exposure time ranged between 30 
and 80 s due to changing sky transparency conditions dur¬ 
ing the nights. We used iraf to reduce the individual im¬ 
ages in a standard way, and to perform aperture photom¬ 
etry on the target and two comparison stars. Differential 
light curves of the target relative to the main comparison 
star (Cl, USNO-A2.0 0900-00554015) were then computed. 
We confirmed the stability of Ci by differential photome¬ 
try relative to the second comparison star. In converting the 
differential 7-band magnitudes of HS 0220-1-0603 to the stan¬ 
dard system, we calculated the 7-band magnitude of Ci from 
its Sloan Digital Sky Survey (SDSS) r, i, and 2 magnitude 
measurements (Data Release 7). We then transformed the 
SDSS magnitudes into 7-band magnitudes by using the two 
formulae given by R. Lupton in 200^ The average value 
for Cl is 7 = 15.565 ± 0.007. 


3.2.2 WHT/ULTRACAM fast photometry 


We also observed HS 0220-1-0603 with ULTRACAM on the 
WHT on La Palma on 2005 August 9 and 10. ULTRACAM 


is a high-speed, triple-beam CCD camera (Dhillon et al. 


20071 capable of acquiring simultaneous images in three dif¬ 


ferent photometric bands. For HS 0220+0603 we used the 
SDSS u', g', and i' filters (see Table 0 . The exposure time 
was fixed to 4.95 s and the dead time between consecutive 


exposures was ~ 24 ms. We reduced the data with the UL¬ 
TRACAM pipeline software. Debiassing, flat-fielding and 
sky subtraction were performed in the standard way. We 
determined the flux of the source and two close comparison 
stars with aperture photometry using a variable aperture 
scaled according to the FWHM of the stellar profile. Note 
that the ULTRACAM observations took place when the sys¬ 
tem’s brightness was starting to increase (see Fig.0, and it 
was about one magnitude brighter than the low state level. 


3.2.3 Updated eelipse ephemeris 


We combined the times of mid-eclipse of the WD in 
HS 0220+0603 in the high state reported by |Rodrfguez-Gil| 
et al. (20071 with new WD mid-eclipse times measured from 


the light curves in the low state presented in this work to 
update the eclipse ephemeris. The resulting linear ephemeris 


ro(HJD) = 2 452 563.57399(12) + 0.149207749(75) x E . 

We used this ephemeris to calculate the orbital phase for all 
the data presented in this work. 


XMM-NEWTON X-RAY AND UV 
OBSERVATIONS 


The XMM-Newton satellite ( Jansen et al.|[^01[ ) pointed to 
HS 0220+0603 for 30.7 ks on 2015 January 13 (observation 
ID 0212480101). The European Photon Imaging Camera 
(EPIC) consists of two MOS cameras ( Turner et al.|[^01[ ) 
and one PN camera (Striider et al. 20011, which were op¬ 


erated in full-frame window imaging mode with the thin 
optical blocking filter applied. We used the Science Analysis 
Software (SAS, v.14.0) to carry out the data reduction and 
obtain the scientific products. The observation was affected 
by episodes of background flaring, therefore a total of 14.2 ks 
had to be neglected in our analysis. We extracted events us¬ 
ing a circular region with a radius of 15 arcsec centred on 
the source position (RA = 02:23:01.6, Dec. = +06:16:49.6, 
J2000), and a circular region with a radius of 30 arcsec cov¬ 
ering a source-free part of the CCDs to extract background 
events. The average 0.2-12 keV source net count rate was 
(8.2 ± 1.0) X10“^ count s“^ for the PN camera. A very lim¬ 
ited number of photons were recorded with the MOS cam¬ 
eras, which did not allow us to perform a spectral analysis. 

We generated the spectrum and the light curve of the 
PN data, as well as the response files, following the stan¬ 
dard analysis thread^ The spectral data were grouped to 
contain a minimum of 10 photons per bin and fit in the 0.3- 
10 keV energy range using XSpec (v.12.8; Arnaud| 19961. 
We assumed a total Galactic Hi column density (Ah) of 
Ah=(3.2 ± 0.1) X 10^° cm~ ^ calculated with the Ah-Av 
relation (Giiver fc 6zel|2009 1 using the Av obtained in Sec- 


tion |3.1| This is consistent with the Ah=4.8 x 10^° cm“^ m 
the direction to the source found by |Kalberla et ar] ( |2005[ ). 
Since HS 0220+0603 is in a low state, we assume that 


http: //www. sdss3. org/drS/algorithms/ 
sdssUBVRITransf orm.php 


® http://xmm.esac.esa.int/sas/current/documentation/ 
threads/ 
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Figure 2. The 0.3-10 keV PN X-ray spectrum (top) and resid¬ 
uals (bottom). The solid line represents the best fit with a one- 
temperature thermal plasma model (PHABS*APEC). 


the bulk of the X-ray emission comes from the M-dwarf com¬ 
panion star. Therefore, we modelled the spectrum with an 
absorbed one-temperature thermal plasma model, which is 
typically used for stellar emission. To this end, the plasma 
thermal emission APEC code (Smith et al. 20011 was ap¬ 
plied. A good fit with 'xf,=Q.7 for 18 degrees of freedom (dof) 
and p-value >0.8 (see Fig.j^ was achieved. The value of the 
obtained temperature is kT = 2.9 keV (~ 34 MK) and the 
emission measure (EM) is 8.98 x 10®^ cm“®. The 0.5-10 keV 
absorbed (unabsorbed) flux is 2.0 x 10“^'* erg cm“^ s“^ 
(2.1 X 10“^"* erg cm“^ Assuming a distance of 780 pc 
(see Section 6.21, the corresponding X-ray luminosity is 
14 X 10^® erg s“ (see Table for spectral results and X-ray 
luminosities for different energy ranges in order to compare 
with the literature). 

Our results are comparable with those obtained for sin¬ 
gle M-dwarf stars ( Schmitt &: Liefke|2004 1, although in our 
case the luminosity is higher (~1 order of magnitude). For 
example, LMC 335 and 2XMM J043527.2-144301 have tem¬ 
peratures of kT = 1.22 keV and kT = 1.22 keV, and 
X-ray luminosities of 1.7 x 10^® erg s“^ (0. 1-7 keV) and 
5 X 10®^ erg s“^ (0.2-2 keV), respectively (Tsang et al. 


2012 Gupta, Galeazzi & Williams 2011|). When placing 


HS 0220-1-0603 in the Lx vs rotation period plot by 
|lato et ^ ( |2003| ) (also see|Cook, Williams &: Berger|| 


Pizzo-| 


20141, 


we see that the luminosity is also higher than rapidly rotat¬ 
ing active stars. We therefore cannot discard the presence of 
X-ray emission due to very low-level accretion. 

Simultaneous to the X-ray measurements, observations 
with the optical/UV monitor telescope (OM; [Mason et al.| 


2001) were carried out in the imaging mode. The near¬ 


ultraviolet uvwl (245-320 nm) and uvm2 (205-245 nm) 
filters were used during ~4.4 ks each. The data were 
reduced using the omichain task. The resulting unred¬ 
dened fluxes for the uvwl and uvm2 bands were (0.24 ± 
erg cm~® s"’^ A \ and (0.48 ± 0.02) x 
respectively (see bottom panel of 

Fig.^. 


0.01) X 10"^® erg cm’® s"^ A ^ 
10"^^ erg cm"^ A 


Table 2. Best 
PHABS*APEC. 

fit to the 0.3—10 

keV PN 

spectrum using 


Parameters model 


Ah 

3.2 xl0®0 cm-2 



kT 

2.9 keV 




8.98 x 10®® cm“® 



X® (dof) 

0.7 (18) 



Energy band 

-^X ?abs 

-^X mnabs 

Lx® 


(10-14 gj.g cm“ 

■® s-l) 

(10®® erg s-i) 

0.5-10 keV 

2.0 

2.1 

14 

0.1-7 keV 

2.03 

2.4 

16 

0.5-10 keV 

1.25 

1.62 

10 


^ Emission measure EM= J* nenfjdV = 4:7vD^NormxF^c- 
To calculate the EM and the luminosities D = 780 pc was as- 
sumed for the distance. 


5 SPECTRAL CLASSIFICATION 


5.1 Spectral energy distribution 

The optical spectrum of HS 0220-1-0603 in the low state is 
dominated by the two stellar components. In Fig. I^the con¬ 
tributions of a WD and a mid-M dwarf are apparent in the 
blue and the red, respectively. The spectral modelling pre¬ 
sented below does not require the presence of an accretion 
disc. In addition, the U-shaped WD eclipse seen in the light 
curves (see Figs.[^and 101 points to a bare WD or, at most, 
a WD surrounded by a cold, remnant disc. The emission 
lines do not reveal a disc either, unlike DW UMa in the low 
state ( jDhillon, Jones fc Marsh||1994| . 

Narrow emission lines of hydrogen are superimposed on 
the WD absorption troughs. For comparison, the intrinsic 
FWHM of Ha in the low state is 340 km s“^ at orbital phase 
0 (defined as inferior conjunction of the secondary star), 
whilst the Ha emission line has a FWHM of 1390 km s“^ in 
the high state ( Rodrfguez-Gil et al.|20()7 l. The Can triplet 
(A8498, A8542, and A8662) in emission is apparent around 
orbital phase 0.5. Some Hel A5876 emission is also present. 
We will later show that these narrow emission lines origi¬ 
nate on the donor side of the binary system, and are due to 
irradiation from a hot WD and/or chromospheric emission. 
The He ll A4686 and Bowen blend emissions observed in the 
high state are absent. 

In determining the spectral type of the two stars we 
used an average spectrum of HS 0220+0603 computed from 
spectra corresponding to a limited range of orbital phases 
around the primary eclipse. This minimises the contribution 
from the side of the secondary star that faces the WD and 
yet includes the spectrum of the WD. In particular, we only 
used VLT/FORS2 spectra obtained at orbital phases 0.15 
and 0.20, WHT/ISIS red spectra at 0.92, 0.93, and 0.11, 
and at 0.82, 0.92, 0.93, 0.11, and 0.21 for the WHT/ISIS 
blue spectra. Note that we use this average spectrum in 
the subsequent analyses unless otherwise stated. Also note 
that this approach will eliminate any information on the 
spectral type variation with orbital phase, but it is necessary 
as the signal-to-noise ratio of the individual spectra is far 
from being optimal for this purpose. 

The spectral contribution of the WD is very blue in 
colour, and contains both Balmer and Hel lines, making 
this the first definite detection of a DAB WD in a CV. The 
composite spectrum superimposed on the average spectrum 
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Figure 3. Average VLT/FORS and WHT/ISIS spectrum of 
HS 0220+0603. Only spectra taken shortly before or after pri¬ 
mary eclipse were used in order include the WD spectrum while 
at the same time minimising the effect of irradiation from the 
WD on the M-dwarf companion. The data (grey) are fitted with 
WD (blue)+M-dwarf (red) templates. The solid black line is the 
best composite fit (see text for details). 


of HS 0220+0603 presented in Fig. [^contains our DAB WD 
model with temperature and surface gravity of 30 000 K and 
logg = 8.35, respectively, and a dwarf M5.5 secondary star 
(see Section 5.21. 

The de-reddened UV fluxes (see Section]^ are plotted 
together with our WD model and the average WHT/ISIS 
blue spectrum in the bottom panel of Fig. The UV ob¬ 
servations fit well into the WD model with TeS = 30 000 
K. For comparison, we show in Fig. the SDSS spectrum 
of a (single) DAB WD with similar atmospheric param¬ 
eters (SDSS J0847 42.22+000647.6, Teg = 30545 ± 241 K, 
logg = 8.49±0.05; Kleinman et al.|2004 | along with the av¬ 
erage spectrum of the WD in HS 0220+0603 and the DAB 
model. The SDSS WD and the HS 0220+0603 spectra show 
a good match. Given the uncertainties in the flux calibra¬ 
tion of the average WD spectrum and the UV fluxes, and 
considering that the UV and the optical observations were 
not simultaneous, it is difficult to establish a formal error on 
the WD temperature, but we conservatively estimate that 
the quoted value is accurate to within 5000 K. 


5.2 Secondary star molecular-band indices 

In this section we present the determination of the spec¬ 
tral type of the secondary star. Boeshaar (19761 intro¬ 
duced a spectral classification scheme for M dwarfs based 
on molecular-band ratios, which has also been used for CVs 


(Wade & Horne 19881. The work of Reid, Hawley Sz Gizis 
1 1995P used a similar set of spectral indices to measure the 
depth of selected molecular band heads. These authors im- 


Figure 4. Top panel: Average WD spectrum around eclipse 
(grey) with our DAB WD model (red; = 30 000 

K, logg = 8.35) and spectrum of the field DAB WD 

SDSS J084742.22+000647.6 (blue; = 30545 K, logg = 8.49) 
superimposed. Note that the contribution of the secondary star 
has not been subtracted, so there is a flux excess increasing with 
wavelength. Bottom panel: Average WD spectrum around eclipse 
(grey) with our DAB WD model (blue) and the XMM UV flux 
measurements after correcting for interstellar reddening (circles). 


proved on the M-dwarf spectral sequence presented by|Kirk-] 


Patrick, Henry Sz. McCarthy (19911 for classification in the 
red (6300 —9000 A). These indices have been extensively cal¬ 


ibrated against spectral type by many authors (see Lepine, 
Rich fc Shara|[2003 and references therein). 


An assessment of the spectral type of the secondary star 
in HS 0220+0603 can therefore be done by measuring a set 
of these spectral indices from the average WHT/ISIS red 
spectrum and comparing them with the indices observed for 
calibrated M dwarfs. We followed the strategy described by 


Lepine et al. (2003 \ for mid to late M-dwarfs, focusing on the 


VOl, V02, and Ti07 indices (see their table 2 for details). 
We did not use the Ti06 index as the continuum is defined 
very close to the O 2 telluric absorption complex at ~ 7650 A, 
and residuals from the telluric correction remain in that re¬ 
gion. Prior to this analysis we first removed the contribution 
of the WD. We present the flux ratios for HS 0220+0603 and 
SDSS M-dwarfs ranging from M2 to M8 in Table |Lepine| 


et al. (20031 present linear functions to obtain the spectral 


type of M dwarfs from the molecular-band indices. Using 
their equations (4), (6), and (7) we obtain a spectral type 
for the secondary star in HS 0220+0603 of M5.5, M5.6, and 
M5.4 for the VOl, V02, and Ti07 indices, respectively. The 
indices obtained for the SDSS M-dwarf templates (Table 
also point to a M5-6 spectral type. The linear calibrations 
of Lepine et al. (20031 can provide spectral types accurate 
to half a spectral type. We can therefore conclude that the 
M-dwarf secondary star in HS 0220+0603 shows a low-state 
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Table 3. Molecular band-head indices. 


Star 

VOl 

V02 

Ti07 

HS 0220-1-0603 

0.875 

0.649 

0.756 

M2 V 

0.947 

0.899 

0.956 

M3 V 

0.920 

0.821 

0.914 

M4 V 

0.904 

0.758 

0.876 

M5 V 

0.874 

0.654 

0.791 

M6 V 

0.851 

0.601 

0.762 

M7 V 

0.829 

0.537 

0.703 

M8 V 

0.799 

0.443 

0.626 



HJD - 2450000 


spectral type of M5.5±0.5 when observed close to its inferior 
conjunction. 


Figure 5. Modelling of the WHT/ULTRACAM WD eclipse in 
the g'-band. The data were fitted with the LCURVE code. The solid 
line represents the best fit. 


6 BINARY PARAMETERS 


6.1 White dwarf eclipse modelling 

In order to obtain constraints on the binary parameters we 
fitted light curve models to the ULTRACAM g'-band eclipse 
data. The g'-band light curve has the best quality and the 
deepest eclipse of the three recorded by ULTRACAM (see 
Fig. §. Of special interest is the mass ratio q = M 2 /Mi 
(where Mi and M 2 are the masses of the WD and secondary 
star, respectively) as a function of the orbital inclination (i), 
that can be tested against an analogous function derived 
from the joint modelling of the emission- and absorption¬ 
line radial velocity curves and the ellipsoidal modulation 
plus eclipse curves (see Section [R^ . 

For the modelling of the ULTRACAM light curve we 
made use of the lcurve code developed by one of us 
(TRM; see Copperwheat et al.||2010| for a description, and 
e.g. IPyrzas et al.||2009| for its application to light curves 
of WD/M-dwarf binaries). The LCURVE code computes a 
model from input parameters provided by the user. Model 
light curves are then fitted to the data using Levenberg- 
Marquardt minimisation. Both fixed and free parameters 
can be set by the user, which gives full flexibility to the 
code operation. 

The results of the spectral modelling presented in Sec¬ 
tion |5.1| allowed us to fix the temperature of the WD to 
Ti = 30 000 K. We also fixed the limb darkening coefficients 


to 4-parameter values as listed in Gianninas et al. (20131 


for Ti = 30 000 K and logg = 8.5. The mass ratio q, the 
scaled WD radius ri = Ri/a (with a the binary separa¬ 
tion) and the effective temperature of the secondary star 
T 2 were left as free parameters. For the scaled secondary 
star radius r 2 = R 2 /CL, the code calculates the proper non- 
spherical shape of the secondary star’s Roche lobe for each 
value of q. We performed a grid search in q, with q ranging 
from 0.25 to 0.60 in steps of 0.05. The results are presented 
in Table m 


6.2 Simultaneous radial velocity and ellipsoidal 
modulation modelling 

The NOT 7-band light curves were phase folded according 
to the eclipse ephemeris given in Section |3.2.3| and aver¬ 
aged into 43 orbital phase bins. To interpret the 7-band light 
curve and the absorption- and emission-line radial velocity 


Table 4. ULTRACAM WD eclipse modelling parameters. 


Mass ratio 
g (= M2/Ml) 

Inclination 
i (deg) 

Scaled WD radius 
ri (= Ri/a) 

0.260 

88.71(15) 

0.0101(5) 

0.300 

85.48(04) 

0.0095(5) 

0.350 

83.51(03) 

0.0091(4) 

0.400 

82.11(03) 

0.0087(5) 

0.450 

80.99(02) 

0.0083(4) 

0.500 

80.05(02) 

0.0080(4) 

0.550 

79.24(02) 

0.0078(4) 

0.600 

78.52(02) 

0.0075(4) 


curves we used the xrbcurve model described in IShahbazI 


et al. (20031, which has successfully been used to model the 


light curves and radial velocity curves of neutron star and 
black hole X-ray binaries ( Shahbaz et al.|2004 L Briefly, the 
model includes a Roche-lobe filling secondary star, the ef¬ 
fects of UV/X-ray heating of the secondary star by a source 
of high energy photons from a compact object, an accretion 
disc and mutual eclipses of the disc and star. In the case 
of HS 0220-1-0603 in the low state, where the eclipse of the 
WD by the secondary star is clearly seen, we simulate the 
presence of a bare WD by assuming that it can be described 
by a blackbody cylindrical accretion disc (with equal radius 
and height). 

The geometry of the binary system is given by the or¬ 
bital inclination i, the binary mass ratio q, and the Roche- 
lobe filling factor /rj of the secondary star. The light arising 
from the secondary star depends on its mean effective tem¬ 
perature (T 2 ), the gravity darkening exponent /3 and the 
X-ray/UV albedo W. The additional light due to UV/X-ray 
heating is given by the unabsorbed UV/X-ray flux {Fx,o)- 
The light from the WD is given by its scaled radius and ef¬ 
fective temperature. The distance to the source in parsecs 
(dpc), the orbital period P and the radial velocity amplitude 
of the secondary star {K 2 ) set the scale of the system. 


We use phoenix model-atmosphere fluxes (Hauschildt, 


Allard & Baron 19991 to determine the intensity distribu¬ 


tion on the secondary star and blackbody fluxes for the light 
arising from the WD. For both cases we use a quadratic limb- 
darkening law with coefficients taken from Claret (19981, to 
correct the intensity. We assume that the secondary star is 
in synchronous rotation and completely hlls its Roche lobe. 
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Table 5. List of the key model variables and what they represent. 


Variable 

Definition 

Ml 

Mass of the white dwarf 

M 2 

Mass of the secondary star 

Q 

Binary mass ratio defined as M 2 /Mi 

i 

Orbital inclination 

p 

Orbital period 

a 

Binary separation 

Fx,o 

Unabsorbed heating flux 

dpc 

Distance in pc 

K 2 

Secondary star’s radial velocity amplitude 

T2 

Secondary star’s mean effective temperature 

Ri 

White dwarf radius 

Ti 

White dwarf effective temperature 

0 

Gravity darkening exponent 

W 

Albedo of the secondary star 

E(B - V) 

Colour excess 


Since the late-type secondary star should have a convec¬ 
tive envelope, we fix the gravity darkening exponent to 0.08 
I Lucy|[l967 l. Its albedo [W) is fixed at 0.40 (de Jong et al. 


1996). 


The irradiating flux can strongly alter the distribution 
of temperature across the face of the secondary star and 
thus the observed optical light curves and radial velocities. 
For the light curves, the heating of the secondary star is 
computed in the same way as described in [Shahbaz et al.| 
(2003), who calculate the increase in unperturbed local ef¬ 
fective temperature due to the irradiating external source. 
Substantial heating of the secondary star also shifts the ef¬ 
fective centre of the secondary, weighted by the strength of 
the absorption and emission lines, from the centre of mass 
of the star. This results in a significant distortion of the 
radial velocity curve leading to spuriously high or low K- 
amplitudes using absorption or emission lines, respectively, 
arising from the secondary star (Phillips, Shahbaz & Pod- 
siadlowski 1999 Shahbaz et al. 2000). To model the radial 
velocity curves we determine the fraction of the surface of 
the star that contributes to the absorption-line or emission¬ 
line radial velocity. These are described by factors Fay and 
Fev, respectively, and represent the fraction of the exter¬ 
nal radiation flux that exceeds the unperturbed flux, e.g. 
Fav= 1-0 assumes all the area of the star visible by the 
heating source is irradiated and does not contribute to the 
absorption-line radial velocity, whereas for Fev=1-10 only 
surface elements on the star where the external radiation 
flux is greater than 10 per cent of the unperturbed flux 
contribute to the emission-line radial velocity (Billington, 
Marsh & Dhillon 1996 Phillips et al. 1999). The line flux of 
each element on the star is calculated using the equivalent 
width and continuum flux value, and the equivalent width 
of elements that do not contribute are set to zero. For ease 
of reference, the key variables are listed in Table|^ 

In determining the binary parameters we simultane¬ 
ously fit the photometric 7-band light curve (43 data points) 
and the absorption-line (Nal doublet) and emission-line 
(Call triplet) radial velocity curves (11 and 7 data points, 
respectively) with our model to represent the photometric 
and radial velocity variations. We used the differential evo¬ 
lution algorithm described in |Shahbaz et al.] ( |2003[ ) to fit the 
data, which is robust and simple. 

The model parameters that determine the shape and 


0.8 



0.2 


0 _I_I_^_I_I_I_ 

0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 

M, (M„) 


Figure 6. Mi—M 2 plane with the best fit solutions. The solid 
red circles show the 99 per cent confidence level solutions. The 
dashed line corresponds to a mass ratio of ij = 0.54. The bars at 
the lower left represent the typical error. 


amplitude of the optical light curves, absorption-line and 
emission-line radial velocity curves are i, q, T 2 , Fx,q, dpc, 
W, Ri, Ti, K 2 , Fay, and Fev- For the fitting procedure 
there are three extra parameters: the phase shift for the 
optical light curve 5(j>i^c, and the phase and radial-velocity 
shift for the absorption-line radial velocity curve Ji/irv and 
7rv, respectively (the absorption-line radial velocity curve 
allows us to define the phase 0.0, i.e. inferior conjunction of 
the secondary star, and the systemic velocity, hence we use 
the same values for the fit to the emission-line radial velocity 
curve). 

We fixed T 2 = 3000 K for a mid-M dwarf (see e.g. 
Leggett et al. [ 19 ^ , a WD blackbody temperature of 


Ti = 30 000 K, E{B — V) = 0.047 with the constraint 
Ml < 1.4 Mq and then performed a grid search in i. Given 
that there are three different data sets with different num¬ 
bers of data points, to optimise the fitting procedure we 
assigned relative weights to them. After our initial search 
of the parameter space, which resulted in a good solution, 
we scaled the uncertainties on each data set (i.e. the light 
curve and the radial velocity curves) so that the total re¬ 
duced of tho fit was ~ 1 for each data set separately. 
The fitting procedure was then repeated to produce the fi¬ 
nal set of parameters. We have checked that the results do 
not significantly change upon varying T 2 by a few hundred 
degrees. 

By combining the best fit model parameters q, i and K 2 
with the mass function equation, we can determine Mi and 
M 2 . In Fig.|^we show the best fit solutions in the {Mi — M 2 ) 
plane, where the most probable value for the binary masses 
are Mi = 0.87 Mq and M 2 = 0.47 Mq. The red circles in the 
figure show the 99 per cent confidence level solutions. Table[^ 
shows the mean and rms uncertainties in the parameters 
values, where we have only used solutions with the 99 per 
cent confidence level. As one can see the (Mi — M 2 ) solution 
allows us to constrain g; the dashed line in Fig. shows the 
best fit with q = 0.54 ± 0.02, which when combined with 
the (i — q) solutions obtained from the ULTRACAM eclipse 
fitting (see Section [6T| | , gives i ~ 79 degrees. In Fig. [To|we 
show the light curve and radial velocity curves for the best 
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Table 6. Results obtained from the model fits to the /-band 
optical light curve, Nal absorption-line and Call emission-line 
radial velocity curves of HS 0220-1-0603 for T 2 = 3000 K, and Ti = 
30000K. 


Parameter 

Value 

g 

0.54 + 0.03 

K 2 (km s“l) 

284 ± 11 

dpc 

740 ± 26 

logFx.o (erg cm~^ 8“^ ) 

-11.7 + 0.05 

7AV (km s-i) 

1.8 + 2.3 

Pav 

1.06 + 0.06 

Fev 

1.03 + 0.02 

Ml (Mq) 

0.87 + 0.09 

M 2 (Mq) 

0.47 + 0.05 



fit model parameters; M\ = 0.87 Mq and M 2 = 0.47 M©, 
i = 79°, Ta = 3000K, and Ti = 30 000 K. 

We also present a model prediction of the Ha emission¬ 
line radial velocity curve (not included in the fitting pro¬ 
cess) assuming that all of the inner face of the secondary 
star is irradiated. The model nicely reproduces the observed 
Ha radial velocity curve, but note that the Ha emission line 
is also seen around zero phase, indicating that it is not en¬ 
tirely driven by irradiation. We suggest that chromospheric 
Ha emission is also occurring in the secondary star (also see 
Mason et al.|20d8 for the case of VV Pup in the low state). 


7 DISCUSSION 
7.1 The white dwarf 

The spectral decomposition of the optical spectrum of 
HS 0220-1-0603 in the low state revealed a DAB WD with 
Tefi ~ 30 000 K and logg ~ 8.35. On the other hand, our 


Figure 7. Comparison of the HS 0220+0603 average H/3 absorp¬ 
tion profile (top grey spectrum) with the H/J absorption profiles 
of three magnetic WDs (showing the characteristic Zeeman split¬ 
ting) with magnetic fields of 8.4, 4.1 and 1.0 MG (from top to 
bottom; see text for details). 


Kulebi et al(][2009t , and SPSS J091005.44+ 081512.2 (Tes = 
25000K, Bi = l.OMG, Kulebi et al.||2009 |. We can there¬ 
fore rule out a WD magnetic field in excess of 8 MG in 
HS 0220+0603, but Bi ^ 4.1 MG would pass undetected 
due to the poor signal-to-noise ratio. Blue spectra of better 
quality are needed to confirm any Zeeman splitting in the 
low state. 

Observations of accreting WDs in GVs indicate that 
they are significantly more massive than isolated WDs. The 
statistical studies of the GV population by [Smith fc Dhillon| 
119981 and Knigge (20061 produced mean WD masses of 


0.80 ± 0.22 Mq (above the period gap, 0.69 ± 0.13 Mq be- 


dynamical solution of the binary system provided a WD 

low) and 0.75 + 0.05 Mq, respectively. Also, 

Savoury et al. 

with mass Mi = 0.87 ± 0.09 Mq. Cubic interpolation of the 

(20111 and Zorotovic, Schreiber & Gansicke 

(20111 report 


q — ri data from the ULTRACAM eclipse modelling (see 
Table 1^ yields a scaled WD radius Ri/a = 0.0078 ± 0.0004 
for a mass ratio q = 0.54, which translates into a WD radius 
Ri = 0.0103 ± 0.0007 Rq. The WD surface gravity derived 
from the measured mass and radius is logg = 8.36, in good 
agreement with the logg ~ 8.35 value found from spectral 
modelling (Section |5.1[ ). When the WD radius is used in the 
WD spectral modelling a distance of ~ 763 pc is derived, 
which is consistent with the results presented in Section 
(see Table [^. 

Circular polarimetry observations of SW Sex stars show 


more massive WDs in GVs. The WD mass for HS 0220+0603 
fits well within these results. In addition, the WD mass of 
the eclipsing dwarf nova IP Peg {P — 3.80 h) is measured to 
be 1.16 ± 0.02 Mq ( [Copperwheat et al.|2010| . On the other 
hand, the mass distribution of isolated WDs peaks at < 0.6 
Mq (e.g. Koester, Schulz & Weidemann 1979 Bergeron, Saf- 
fer & Liebert |1992| K epler et al. [2007 Falcon et al. |2010| 
Tremblay et al. 20131. 


that some systems may contain magnetic WDs (Rodriguez- 


Gil et al. 2001 2002 20091. In order to check for any Zeeman 


splitting in the WD spectrum of HS 0220+0603 we compare 
the average H/3 absorption profile of HS 0220+0603 with 
three magnetic WDs (Fig. [t]): SDSS J151130.17+422023.00 
(Teff = 30882K, Bi = 8.4MG, [Kepler et al.||2013|), SDSS 
J154305.67+343223.6 (Tett = 25 000 K, Bi = 4.1 MG, 


The distribution of WD effective temperature as a func¬ 
tion of the orbital period may serve as a probe of CV evo¬ 
lution. Changes in the mass transfer rate reflect in the WD 
effective temperature, which is determined by compressional 
heating of the WD by the accreted matter. Therefore, Tefi 
can be a good tracer of the medium-term (10® — 10® yr) 
mass transfer rate ITownsley fc Bildsten|[^03[ Townsley & 


Gansicke 20091. The WD effective temperatures predicted 
by the evolutionary sequences of [Knigge et al.| ( |2011[ ) are 
significantly cooler than the data below the period gap, and 
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Figure 8. Period-mass data of CV donors (data from |Knigge| 
|2006| |. All mass estimates come from superhump period excesses. 
The SW Sex systems are marked with grey squares. The large 
grey square corresponds to our results for HS 0220+0603 while 
IP Peg is plotted as a large grey triangle. The dark grey solid 
line shows the best-fit evolutionary model track oflKnigge et al.l 

i |2MT| i. 


show a large disagreement in the 3-4 h orbital period range, 
where the scarce measurements point to much hotter WDs 
than expected (see fig. 5 of Townsley fc Gansicke|[2009 and 
fig. 16 of Knigge et al. 12011 1 . However, the relationship be¬ 
tween the time-averaged mass transfer rate and the WD 
temperature has to be used with caution. As [Townsley fc| 


Gansicke (20091 point out, the mass transfer rate for a given 


Tefi strongly depends on the WD mass, which again high¬ 
lights the need for independent mass measurements. In ad¬ 
dition, a dominant WD in the blue/UV is also a requirement 
for a reliable determination. Hence, studies of the WDs 
in CVs in the low state like the one presented here are im¬ 
portant to test the current CV evolutionary sequences in 
the 3-4 h regime. In this regard, Knigge et al. (20111 sug¬ 
gest that the large scatter of the WD temperatures observed 
in this orbital period interval may be the result of the tem¬ 
peratures not tracing the secular mass transfer rates. 

The measured temperature and mass of the WD in 
HS 0220-1-0603 result in a time-averaged mass transfer 
rate of (M) = (1.1 ± 1.2) x 10“® Mq yr“^. This value 
places HS 0220-1-0603 closer to the evolutionary sequences of 
Knigge et al. ( |2011 1 than MV Lyr, TT Ari and DW UMa, 


but at the same time increases the scatter in M observed 


close to the upper edge of the period gap (fig. 5 of Towns- 
ley fc Gansicke||2009 1. In contrast, IP Peg (also in the 3-4 


h orbital period regime) may have a time-averaged mass 
transfer rate as low as (M) < 5 x lO”'^^ M© yr“^ as de¬ 


rived from its WD temperature (10 000 — 15 000 K; Copper- 
jwheat et al.|[MTo| , almost 30 times lesser than the value 
inferred for HS 0220-1-0603 from its WD temperature. These 
authors also point out the possibility that the WD temper¬ 
ature might not be a good tracer of the medium-term mass 
transfer rate. The time-averaged mass transfer rate of IP 


Figure 9. Same as Fig. but for the donor mass—spectral type 
data. HS 0220-1-0603 and IP Peg largely deviate from the evolu¬ 
tionary track ofjKnigge et al.|l|2011^. 


Peg is also inconsistent with the results of the population 
synthesis study by Howell et al. (20011. It is therefore clear 
that more direct measurements of masses and temperatures 
of WDs in the 3 — 4 h orbital period regime are needed. 


7.2 The M-type companion 


We have shown that HS 0220-1-0603 contains a M-dwarf 
donor star with mass M 2 = 0.47 ± 0.05 M© (Table |^. 
Note that the fits presented in section |6.2| never produced 
a secondary star mass smaller than 0.3 Mq. We have ap¬ 
proximated the radius of the donor star by its Roche-lobe 
volume radius | Eggleton||1983 1. The donor radius can then 
be derived from its scaled value R 2 /a by using the binary 
parameters of HS 0220-1-0603. The resulting Roche-lobe vol¬ 
ume radius of the M-dwarf is then R 2 — 0.43 ± 0.03 R©. 
The distance to the system when this radius is taken into 
account in the M-dwarf spectral modelling is found to range 
between 740 and 940 pc for the spectral types M5 and M6, 
respectively. This is consistent with both the distances from 
the WD spectral modelling and the dynamical modelling 
presented in Section [6.2| 

However, our derived mass for the secondary star runs 
counter to expectations for an orbital period of 3.58 h. The 


semi-empirical CV evolution track of Knigge et al. (20111 
points to an M-dwarf mass of 0.25 M©. In Fig.j^we show the 
P-M 2 diagram for CVs along with the evolutionary track of 


Knigge et al. (20111. Note that all mass values are not direct 


measurements but come from superhump period excesses, 
as emphasised in Section The SW Sex systems between 3 
and 4 hours have been plotted as grey filled squares. We have 
also plotted the data for the dwarf nova IP Peg (large grey 
triangle), with an orbital period of 3.80 h and a secondary 
star mass of 0.55±0.02M©. The secondary star in IP Peg is 
significantly more massive than predicted by the evolution¬ 
ary track of jKnigge et al.j ( |2011| ). In addition, the measured 
mass of the secondary star in HS 0220+0603 marginally fits 
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in the predicted orbital period-secondary mass distribution 
for CVs of Howell et al. ( 2001| (cf. their fig. 5), while the 
secondary mass in IP Peg is clearly off. 

With measured values of 0.55 Mq and 0.47 R© |Cop-| 


perwheat et al. (20101 conclude that the secondary star in 


IP Peg is undersized for its mass, suggesting it is in thermal 
equilibrium. The same might be happening in HS 0220+0603 
with a secondary star mass of 0.47 M© and a derived sec¬ 
ondary Roche-lobe volume radius in HS 0220+0603 of 0.43 
R©. Taken at face value, this indicates that the donor star 
in HS 0220+0603 may not be bloated at all. This is not com¬ 
mon among CVs, in which the secondary stars are believed 
to be significantly oversized for their mass as the result of 
either nuclear evolution or being out of thermal equilibrium. 

In Fig. we have plotted the (few) CVs with determi¬ 
nations of both the secondary star mass and spectral type. 
HS 0220+0603 is marked by a large grey square and IP Peg 
is plotted with a large grey triangle. Our analysis provided 
a secondary star spectral type of M5.5+0.5 when the donor 
star is observed at its inferior conjunction (i.e. phase zero). 
This indicates a photosphere with temperature around 2800 
K (see e.g. Bonnefoy et al. 2014| and references therein). 
Even if HS 0220+0603 followed the semi-empirical sequence 
of Knigge et al. (20111 and had a donor mass of 0.25 M©, the 
spectral type would still be cooler than expected. Fig. I^also 
shows that the secondary star in IP Peg is underluminous for 
its mass. The spectral types measured for HS 0220+0603 and 
IP Peg are also much later than predicted by [Howell et al.| 
(20011 (cf. the effective temperature-orbital period distribu¬ 
tion for CVs in their fig. 7). Again, we stress the fact that 
we measured the M5.5+0.5 spectral type of the secondary 
star in HS 0220+0603 from the average of spectra taken at 
orbital phases around its inferior conjunction, so any infor¬ 
mation on spectral type change along the orbit is lost. An 
illustrative example is the case of the polar CV ST LMi, 
whose secondary star shows a changing temperature with a 
much cooler region at orbital phase 0.8 ( [Howell et al.|2000| ). 
Therefore, more detailed studies of the secondary stars in 
the nova-like variables which populate the 3-4 h orbital pe¬ 
riod range during low states should be conducted in order 
to check how frequent these secondary stars are. 


8 CONCLUSIONS 

We have presented the first dynamical characterisation of 
an eclipsing SW Sex star in the 3-4 h orbital period range 
during a low state (see Table [^. 

We obtained time-resolved optical spectroscopy and 
photometry of HS 0220+0603 during its 2004-2005 low 
state. We refined the eclipse ephemeris by including new 
mid-eclipse times form the photometry. Spectral modelling 
of the secondary star spectra taken around primary eclipse 
revealed a spectral type of M5.5+0.5, suggesting an effective 
temperature of about 2800 K. The analysis of the spectrum 
of HS 0220+0603 also provided the first detection of a DAB 
WD in a CV, with = 30 000 ± 5000 K. By modelling 
the ULTRACAM light curve of the WD eclipse we deter¬ 
mined solutions in the binary mass ratio-inclination angle 
plane. Furthermore, by simultaneously fitting the radial ve¬ 
locity and irradiated ellipsoidal light curves of the secondary 
star, we independently calculate solutions as well as the un¬ 


Table 7. System parameters for HS 0220+0603. R 2 is the volume 
radius of the secondary star’s Roche lobe as defined by [Eggletonj 
( |1983| |. 


Parameter 

Value 

9 

0.54+0.03 

a (R©) 

1.31 ± 0.03 

Ml (M©) 

0.87+0.09 

i?l (R©) 

0.0103 + 0.0007 

Ti (K) 

30 000 + 5 000 

M 2 (M©) 

0.47+0.05 

R 2 (R©) 

0.43 + 0.03 

K 2 (kms“^) 

284 + 11 

i n 

~ 79 

d (pc) 

740 + 26 


perturbed radial velocity amplitude of the secondary star 
{K 2 = 284 ± 11km s“^) and the distance to the source 
(d = 740 + 26 pc). Combining the results of the white-dwarf 
and secondary-star modelling we derived the stellar masses 
to be Ml = 0.87+0.09 M© and M 2 = 0.47+0.05 M© for the 
WD and M dwarf, respectively. As in IP Peg, the M-dwarf in 
HS 0220+0603 seems too cool and undersized for its mass. 
More binary parameters from dynamical measurements in 
nova-like CVs with orbital periods between 3 and 4 hours 
during low states are needed before a critical comparison 
with theoretical predictions can be made. What triggers low 
states in CVs is still far from being fully understood. 
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HS0220 + 0603 — low state 
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Figure 10. From left to right and from top to bottom: NOT /-band ellipsoidal-|-WD eclipse light curve; Call emission triplet radial 
velocity curve; Nal absorption doublet radial velocity curve; Ho emission radial velocity curve. Solid lines are the best fits with the 
exception of Ho (model-predicted curve). The plots under each panel are the corresponding residuals. The orbital cycle has been 
repeated for clarity. 


isation for Astronomical Research in the Southern Hemi¬ 
sphere, Chile, under program 074.D-0657, and on observa¬ 
tions obtained with XMM-Newton, an ESA science mission 
with instruments and contributions directly funded by ESA 
Member States and NASA, and on observations made with 
the William Herschel Telescope operated on the island of La 
Palma by the Isaac Newton Group in the Spanish Observa- 
torio del Roque de los Muchachos, and on observations made 
with the Nordic Optical Telescope operated by the Nordic 
Optical Telescope Scientific Association at the Observatorio 
del Roque de los Muchachos, La Palma, Spain, and on ob¬ 
servations made with the IAC80 telescope operated on the 
island of Tenerife by the Institute de Astrofisica de Canarias 
in the Spanish Observatorio del Teide. The William Herschel 
Telescope data were obtained as part of the 2004 Interna¬ 
tional Time Programme of the night-time telescopes at the 
European Northern Observatory. 
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